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Abstract 

I It is shown a relation between the Fick-Jacobs equation and Schro- 

dinger equation. The Fick-Jacobs equation describes the diffusion in 
a channel with a shape of a surface of revolution. For the case of con- 
stant diffusion coefficient the following results are obtained: the case 
of the conical channel is mapped to quantum free particle; if the shape 
of channel is a throat is obtained a relation with quantum particle in a 
constant potential; also the sinusoidal channel is related with quantum 
particle in a constant potential; in addition, when the channel cross 
section varies as a quadratic exponential, is shown its equivalence with 
the quantum harmonic oscillator. The general case of variable diffu- 
sion coefficient is also considered. In this case a change of variable 
is proposed, and if it is invertible, the Fick-Jacobs equation is equiv- 
alent to the Schrodinger equation of a particle in an effective potential. 
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1 Introduction 

Recently, mathematical techniques used to study certain physical systems 
have been used to address another physical systems that are completely 
different. For example, with the so called AdS/CFT correspondence, us- 
ing theories of gravity it has been able to find results in field theories pQ. 
Moreover, with this duality it has been possible to relate condensed mat- 
ter systems with gravitational systems [21 E]. In addition, Unruh found an 
analogy between the theory of fluids and some spacetimes like black holes 
[I], other works on this subject can be seen in [5]. Is worth stressing that 
by studying some metamaterials it can be simulated certain spacetimes [5J. 
Also, recently methods of physics haven been used to address some problems 
in other disciplines such as finance and economics. For instance, the Black- 
Scholes equation that plays an important role in finance, can be mapped 
to the Schrodinger equation of a free nonrelativistic particle [TJ. Is worth 
noticing that the Hamiltonian obtained in such case has a supersymmetric 
partner [8]. With such Hamiltonians one can propose a generalization of 
supersymmetric quantum mechanics [9]. 

Now, there are many problems related to the transport of particles through 
channels or other confined system. The most simple way to study them is 
through the Fick diffusion equation [TJ], that can be mapped to Schrodinger 
equation. However in several cases one need to use more general equations to 
describe the problem [TUjlT]. Particularly, when the channel has the shape 
of a surface of revolution whose cross section has area A(x), the problem re- 
duces to solving an effective diffusion equation in one dimension known as the 
Fick- Jacobs equation [T2J[I3]- This is an effective one dimensional equation 
that approximates the diffusion on higher dimensions for confined systems by 
introducing an effective diffusion coefficient that varies with position D(x): 



where C(x,t) is the concentration of particles that diffuse. Notice that if 
the area and the diffusion coefficient are constant, we obtain the usual Fick 
equation. In standard references of Fick- Jacobs diffusion [TOj, El [13] one of 
their main goals is to find the corresponding relaxation time of the system. 
This is achieved either by analytical and perturbative methods or through 




(1) 



2 



numerical simulations. 

In this work is shown a relation between the Fick- Jacobs diffusion equa- 
tion and the Schrodinger equation. It is shown that such relation can be 
realized in several cases. When the diffusion coefficient is constant the fol- 
lowing results are obtained: the case of the conical channel is mapped to 
quantum free particle; if the shape of channel is a throat, we have a relation 
with quantum particle in a constant potential; also the sinusoidal channel 
is related with quantum particle in a constant potential; in addition when 
the channel cross section varies as a quadratic exponential, we obtain a re- 
lation with the quantum harmonic oscillator. For the general case when the 
diffusion coefficient is not a constant but a function of the position D(x), it 
is found an equivalence between Fick- Jacobs and Schrodinger equations, if 
a certain change of variable is invertible. The result is that the Fick-Jacobs 
equation with variable diffusion coefficient can be mapped to the Schrodinger 
equation for a nonrelativistic particle in an effective potential. 

This paper is organized as follows: In section 2 is presented a review of 
quantum mechanics with a generalized momentum operator. In section 3 
is studied the Fick-Jacobs equation for the conical channel. In section 4 is 
considered the general case with constant diffusion coefficient and are also 
presented three cases with exact solution. In section 5 is analyzed the general 
case. Finally in section 6 the summary is given. 

2 Generalized quantum mechanical momen- 
tum operator 

In this section we introduce a generalized momentum operator with which 
one can define more than one single Hamiltonian; in particular, a pair of su- 
persymmetric Hamiltonians and two non-Hermitian Hamiltonians. With the 
non-Hermitian Hamiltonians we can obtain wave equations that are equiva- 
lent to Schrodinger equations. This section describes the method used later 
to demonstrate that the Fick-Jacobs equation is equivalent to a Schrodinger 
equation. 
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Lets consider the following Hermitian operator 

that is the usual definition of momentum operator in quantum mechanics 
with h = 1. We can define a different non-Hermitian operator by adding the 
derivative of some arbitrary function f(x) to the original momentum giving: 

PM = eftoPe-™=P + i?&L. (3) 
This generalized momentum operator was first studied by P.M. A Dirac [15J. 

With the new momentum operator ([3]), we can build the following two 
Hermitian Hamiltonians: 

'2 r / ic\ 2^ 



* - « 2 ^>= a2 P 2 +iH£i I 



* = (I) J ■ (5) 

where a is a constant. Notice that if we substitute in (jj]) and (E]) the following 
definition 

f(x)= f X W(u)du, (6) 



o 



we obtain 



With these two Hamiltonians one can construct a supersymmetric quantum 
mechanics |16j . Furthermore, with the generalized momentum operator ([3]) 
we can also build two non-Hermitian Hamiltonians: 

^ ^.^ 2 (^- 2 #-S-(I) 2 )' < 9 > 



where (5 is a constant. 



Lets notice that the transformation f(x) — > —f(x), implies that if 3 — > 
H4, which means that once we solve one problem, we will have the solution 
of the other. It should be noticed that this kind of Hamiltonians arise in 
the so-called quantum finance [7], these Hamiltonians allow to generalize the 
supersymmetric quantum mechanics [9]. A generalized quantum mechanics 
with non-Hermitian operators is studied in |17j . 

With the Hamiltonian H% and H4 we can write the following wave equa- 
tions 

= HsM*,t), (ID 
= &«.(*, t). (12) 

where, as usual if)j, with j = 3,4, is the wave function. The Hamiltonians 
(Q and ({TO]) are non-Hermitian, but the wave equations (flT!) and (fl2l) are 
equivalent to a nonrelativistic free particle. To prove this claim first note 
that 

pt_ pt _ e -f(x) p e f(x) e -f(x) p e f(x) _ e -f(x) p2 e f(x) ^ 

p (/) p (/) = e mp e -f(*) e mp e -m = e f{x) P 2 e- f{x \ (13) 

which corresponds to the equations 

i W&'t) = pr 4 ^ 4 ( x ,t) =/ 3V^P 2 e-^V 4 (x,t), 
from where we have 



d (e f Wip 3 (x,tj) . , „ , 

i-± 1 = /3 2 P 2 (e^V 3 (M) 
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both equations represent the wave equation of a free particle with corre- 
sponding wave function e ± ^ipj- The general solutions of these equations are: 



e 

e 



™M*> t) = e~^ pH (e-^Vo4(*)) , (14) 

where i/; 3 (x,0) = i[) 03 (x), ^(x, 0) = ip 04: (x), are the initial conditions. Thus, 
the general solutions of the equations ffTTj) and (IT21 are given by: 

^ 4 (x, t) = (e f ^e~^ 2p2t e- f ^) tfj 0A (x). (15) 

Similarly, one can prove that the wave equations with a potential V(x) 

.(hfai^ = {jpp^ +v{x) } MXit) (16) 

= (^P {f) P {f) + V(x))Mx,t), (17) 
are equivalent to the following Schrodinger equations 

i-* |A^Z = H (e^Ux,t)), (18) 



(19) 



where the Hamiltonian is just 

H = fi 2 P 2 + V(x). (20) 
The general solutions of the equations (TIB"]) and ( fT9l) are respectively: 

e^Vs(M) = e-^(e^Vo 3 (x)), (21) 

e-WMx>t) = e ~ lAt (e- fix) Mx)) , (22) 
therefore, the general solutions of the equations (fl6l) and ( TT71) are given by 

VaOM) = (e-^e-^V^)^03(a:), (23) 

if> 4 (x,t) = (e /(x) e- ii? *e- /(x) )^04(a;), (24) 



where the initial conditions are i/j 3 (x,0) = i[) 03 (x),i/;4 : (x,0) = ip 04: (x). 



In the next section we will use these results to show that the Fick- Jacobs 
equation is equivalent to a Schrodinger equation. 



3 Fick-Jacobs equation for a conical channel 

In this section we show that the Fick-Jacobs equation for the case where 
the diffusion coefficient is constant and the channel has a conical shape, is 
equivalent to a Schrodinger equation of a nonrelativistic free particle. The 
conical tube was first solved in [11] with some particular choice of boundary 
conditions, here we solve the equation in general. 

If the channel is a cone, the cross-sectional area is A(x) = 7r(l + Xx) 2 , 
where A is the slope of the cone, and if D(x) = D Q is a constant, the Fick- 
Jacobs equation ([I]) takes the form 

dC(x,t) f d 2 C(x,t) 2X dC(x,t) 2A 2 \ 
-^ = D »{—* 1 + Ax dx + (IT^ (M) J- (25) 

Moreover, one can show that taking /3 2 = D Q and f(x) = ln(l + Ax) this 
equation is equivalent to the Schrodinger wave equation: 

-?%^ = H 4i C{x,t). (26) 
Indeed, lets notice that making the identifications 

df{x) X d 2 f{x) -A 2 (df{x)\ 2 A 2 



dx "l + Ax' dx 2 '(1 + Ax) 2 ' \ dx ) ~"(1 + Ax) 2 ' 
the Hamiltonian H^, given by equation ({TO]) , can be written as follows 

H, = D (p 2 + 2i^—P 2A ~ ^ 



(27) 



1 + Xx (1 + Aa;) 2 J 

D (-*- + 2^^ - 2X2 \ (28) 
°\ dx 2 I + Xxdx {l + Xx) 2 )' 1 ; 

With expression ( l28l) . we can see that equation ( |26|) is equivalent to (l25j) . 
In section 2, we showed that equation ( 1261) is equivalent to a nonrelativistic 
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free particle. Therefore, the Fick- Jacobs equation with constant diffusion 
coefficient in a conical channel is equivalent to the Schrodinger equation of a 
nonrelativistic free particle. 

Using the results of Section 2, we can find the solution of the equation 
(F22D to be 

C(x, t) = e f{x) e- Dop2 e- f{x) C (x), (29) 

which satisfies the initial condition C(x, 0) = Co (a;). Considering the explicit 
form of f{x) one obtain 

C(x,t) = (l + \x)e- D ° p2 (^y (30) 

Notice that if 

C (x) = l + \x, (31) 

then 

C(x,t) = l + \x. (32) 
Thus, the system does not evolve. 

Here we presented the relation between Schrodinger and Fick- Jacobs 
equations for a special geometry of the channel. However, the form of the 
channel can be different. In the next section we explore different geometries 
for the channel shape when the diffusion coefficient is constant. 

4 Fick-Jacobs equation with constant diffu- 
sion coefficient 

Motivated by the resulting mapping of the previous section, lets consider the 
case of channel with shape of a surface of revolution and when the diffusion 
coefficient is constant. In this case the Fick-Jacobs equation is equivalent to 
a Schrodinger equation and we study three cases that can be solved exactly. 
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To do so lets consider a constant diffusion coefficient D(x) = D = const., 
such that the Fick- Jacobs equation (pQ) takes the form 



dC{x,t) 
dt 



d 2 d\nA(x) d fd\nA(x)\ 2 1 d 2 A(x) 



dx 2 



dx dx 



+ 



dx 



A(x) dx 2 



C(x,t). (33) 



In the same way as in the previous section, we can now find a function f(x) 
and a effective potential V(x), such that, we can construct the corresponding 
wave equation. Identifying 



, 1. , , 11 d 2 A{x) 1 fd\nA{x) 

fix) = -lnA(x), V(x) = — — - - \ ' — — 

Jy ' 2 y h y ' 2A(x) dx 2 4 V dx 



(34) 



the equation (I3"3l can be rewritten as 

_ dC(x,t) =D J P 2 + 2i df(x)_ p + d 2 f(x) _ f df(x) \ 2 
dt \ dx dx 2 \ dx I 



+ V(x) C{x,t).(35) 



Using the results of Section 2 we can see that this equation is equivalent to 



_ dC{x,t) f(x) 



dt 



e f{x) He- f(x) C{x, t), H = D (P 2 + V{ 



x 



(36) 



which in turn is equivalent to 

d(e- f ^C(x,t) 



dt 



He- f ^C(x,t). 



The general solution of this equation is 

e- f{x) C(x,t) = e-^e-^Coix), 



(37) 



(38) 



where the initial condition is C(x, 0) = Cq(x). Therefore, the general solution 
of the Fick- Jacobs equation with constant diffusion coefficient D(x) = Do is 



C(x,t) = (e f We- 6t e- f W) C Q (x), 



namely 



C(x,t) = jA(x)e 



-Ht 



Cp{x) 
A(x) 



(39) 



(40) 



We can see, that if 

C {x) = ^A(xj, (41) 

then 



C(x 1 t) = y /A^j. (42) 
Thus, the system does not evolve. 

Now consider the following cases: 

4.1 Throat like channel 

If the channel where the diffusion takes place is shaped like a throat, the 
cross-sectional area can be taken as an exponential function 

A(x) = e ax+ P, (43) 

where a and /3 are arbitrary constants. In this case the potential defined in 
HD is 

V(x) = ^. (44) 
Therefore, in this case the general solution is 

C(x, t) = e-^e^e-^V^CVx). (45) 

4.2 Sinusoidal channel 

We can also obtain the solution of a channel with a different geometry, which 
has been studied in the literature, we mean the sinusoidal channel [ID], for 
which the cross section is 

A(x) = B(smjx) 2 , (46) 



being B and 7 constants. In this case the general potential in ( 1341) is simply 

V(x) = -7 2 . (47) 

Then the solution is 

C(x, t) = e D ^ H sin( 7 x)e-^ 2t -^-. (48) 

sm(7x) 
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4.3 Gaussian channel 

If the area of the cross section of the channel where the diffusion takes place 
is shaped like a Gaussian, we have 

A(x) = e ax2+bx+c , (49) 

where a, b and c are constants. The potential defined in (j34|) in this case 
turns to be 

V(C) = a 2 ( 2 + a, with ( = x + -. (50) 

a 

Hence, the Hamiltonian takes the well know form 

H = D (P 2 + a 2 C + a) , (51) 

that is the Hamiltonian of an harmonic oscillator, as long as we identify the 
mass m = 1/(2D ) and the frequency u 2 = AD^a 2 . In this case, the general 
solution is 

C(x,t) = e-^e^e-^+^K-^CaiC). (52) 



5 General Fick-Jacobs equation 

In this section we study the general version of the Fick-Jacobs equation, 
we use a change of variable for which this equation can be written as a 
Schrodinger equation in a effective potential. 



The Fick-Jacobs equation ([!]) can be written as 



dC(x,t) d_ 
dt dx 



^ d 2 C(x,t) d 



ln 



D(x) 
A{x) 



9 ( nt ^lnA(x)\ 

Lets introduce the following change of variable 

dz 



y 



JD z) 



Xq = const., 



dC{x,t) 
dx 



(53) 



(54) 
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notice that with this change of variable we have 



d_ 

dx 







'd{x) d v 

Then fl53l) takes the form 



D(x. 



d 2 



d 2 d(JD(x)) d 



dC(y,t) d 2 C(y,t) 



Of 



dx 2 dy 2 



d 



dx 



D(x) 



A(x) 



dy' 



dC(y,t) 
dy 



(55) 



(56) 



If the change of variable (|54p is invertible, each of the terms of this equation 
can be written as a function of the new variable y. In this case we can define 
the next functions: 



1^5 



In 



D(x) 
Mxj~ 



V(y) 

then (1561) can be written as 

2 



[df(y)\ 2 _ d 2 f(y) + 8_ ( d\nA{x) 
\ dy J dy 2 dx \ dx 



(57) 



dC(y,t) 
dt 



_ 2 dmi_ _ &m + (dm) _ v{ ) c{ ty 

dy 2 dy dy dy 2 \ dy J 



Now, with the definitions 



\ dy J 



we find that C(x,t) must satisfy 

dC(y,t) 
dt 



H f C(y,t). 



P = -i^-, (58) 
dy 



(59) 



Moreover, using the results of Section 2, we obtain 

H f = e ny) He~ f{y \ H = P 2 + V(y). 



(60) 
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Therefore, the equation (1591) can be written as 



= efWHe-MC(y,t), (61) 



from where 

d(e-fMC(y,t) 



H(erWC(y,tj), (62) 



dt 

the solution of this equation is 

e-MC{y, t) = e- tit (e-^C {y)) . (63) 

Therefore, if the change of variables (I54j) is invertible, the general solution of 
the Fick- Jacobs equation is 

C(y t t) = (efMe- At e-™)C (y), (64) 
which satisfies the initial condition C(y,0) = C (y). 



The results of this section are valid only if the change of variable (Jo! 
is invertible. There are several examples where this change of variable is 
invertible, for instance D{x) = x n , n = 1, 2, • • • , or D(x) = e ax , a = const. 

In some interesting cases [1~0| [13], general diffusion coefficient is of the 
form 

= ^ (65 ) 



1 + 



dx 



In that case depends only on the shape of the channel if the change of variable 
( 154]) is invertible or not. For instance, consider the case of the conical channel 
with cross section A(x) = n(l + Ax) 2 , the corresponding change of variable 
is 

which is clearly invertible. We will need to take into account the invertibil- 
ity of the change of variable, case by case, for fl65|) it only depends on the 
corresponding cross section area. 
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6 Summary 



We have provided a connection between Fick- Jacobs diffusion equation and 
the Schrodinger equation. The Fick- Jacobs equation is a one dimensional 
effective equation with variable diffusion coefficient that describes the diffu- 
sion in a channel with a shape of a surface of revolution. To do so first we 
introduced a generalized momentum operator with which one can construct 
four different Hamiltonians, two of them are Hermitian and the other two are 
non- Hermit ian. In the case of constant diffusion coefficient we obtained the 
following results: the case of the conical channel was mapped to quantum 
free particle; in the case when the shape of channel is a throat, is obtained 
a relation with quantum particle in a constant potential; also the sinusoidal 
channel was related with quantum particle in a constant potential; in ad- 
dition when the channel cross section varies as a quadratic exponential, we 
obtained a relation with the quantum harmonic oscillator. The general case 
of variable diffusion coefficient was also considered. In this case we proposed 
a change of variable, if it is invertible, the Fick- Jacobs equation is equivalent 
to the Schrodinger equation of a particle in a effective potential. 
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